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O ■ Abstract 

<^> ' Current experiments at the B factories, designed to perform pre- 

""^ ■ cision measurements of matter-antimatter asymmetry in the B meson 

1} , system, have a much broader physics reach especially in the sector 

Q^j of quarkonium spectroscopy. Here we present a minireview on the 

new charmonium-like states observed at the B factories including the 
X(3872) and Y(4260). 



X 

b : 1 Introduction 



More than 30 years after the November Revolution, quarkonium spectroscopy 
seems to be on a great resurgence trail thanks to the discovery of many new 
states at the B factories and elsewhere. It all began with the discovery of the 
enigmatic _D s j(2317) resonance by the BABAR Collaboration in April 2003[1]. 
Soon after CLEO, Belle, CDF and D0 at Fermilab, and BES joined the 
marathon raising the number almost to the double digits. This collection 
of new states has in some cases confirmed quark model and lattice QCD 
calculations while a few others poses serious challenge to our current under- 
standing of the strong interaction. In this report I present an experimental 
survey of new charmonium-like states discovered and studied at the B facto- 
ries; the BABAR experiment [2] at the PEP-II asymmetric energy e + e~ storage 
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rings located at the Stanford Linear Accelerator Center (SLAC) in the US, 
and the Belle experiment [5] at the KEKB energy-asymmetric e + e~ collider 
in Tsukuba, Japan. 

This minireview starts with some brief remarks on unique features of the 
experiments those in conjunction with the high precision data delivered by 
the B factories have made this quarkonium sojourn a reality! The bulk of 
this report is concentrated on the six new states, namely A(3872), X(3940), 
F(3940), F(4260), F(4320), and Z(3930) followed by possible theoretical 
model(s) to explain them. In the last section I summarize my conclusions 
about these states. 

2 Experiments at the Beauty Factories 

Both 5-factory experiments, BABAR and Belle are designed to carry out preci- 
sion CP violation measurements on the beauty meson sector, so it is natural 
they have so much in common. In particular, they have: 

• A precision tracking system composed of a silicon vertex detector, 
whose main role is to locate the secondary decay vertices close to the 
origin, and a central drift chamber that measures momenta and angles 
of charged particles. 

• Excellent charged-particle identification (PID) capability accomplished 
by combining information on the specific ionization (dE/dx) in the two 
tracking devices with that of a dedicated PID system. For the latter, 
BABAR uses an innovative design in the Detector of Internally Reflected 
Cerenkov radiation. Belle makes use of an array of time-of-flight coun- 
ters whose performance is augmented by an Aerogel Cerenkov Counter 
system at higher momenta. 

• An electromagnetic calorimeter, comprising CsI(Tl) crystals, which is 
used to measure the energies and angular positions of photons and 
electrons. 

• Iron flux return of the magnet instrumented with resistive plate cham- 
bers to identify muons and neutral hadrons over a broad kinematic 
range. 

Figure [1] shows the layout of the BABAR and Belle detectors. Although it is 
the overall performance that drives the thing, the most important detector 
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Figure 1: Layout of the BABAR and Belle detectors. 



components for spectroscopic studies are the PID system and the calorime- 
try. To add on to this brilliant detector performance, the B factories are 
accumulating high quality data by leaps and bounds. The results reported 
here are based on data collected by BABAR and Belle although relevant results 
from other experiments will be cross-referenced as appropriate. 

3 Ensemble of Exotica 
3.1 X(3872) - the poster boy 

The X(3872) was first observed in the decay B~ -> X(3872)-FC~, X -> 
J /if) 7T + 7r _ [l] by the Belle Collaboration [5] and subsequently confirmed by 
the BABAR®, CDF [7] and D0[8] Collaborations. The mass of this state is 
m = 3871.2 ±0.5 MeV and the width is T < 2.3 MeV at 90 % confidence level 
(CL) which is consistent with the detector resolution. Study of the 7r + 7r~ 
invariant mass distribution in the decay X(3872) — > J/ip tt + tt~ suggests that 
it may proceed through an intermediate p° resonance. If so, one can expect 
to find its charged isospin partner X(3872)~. BABAR has searched for this 
state in the decays B~ — > J/ip tt^tt Kg and B° — > J/ip tt~7t K + . However, no 
evidence for a charged X(3872) has been found[9] and we set the following up- 
per limits at 90% CL: B(B~ -> X(3872)~K^, X~ -> J /if) n'Tr ) < 11 x lO" 6 
and B{B° -»■ X{3872)~K + , X~ ->■ J/ipTi-jT ) < 5.4 x 10~ 6 . A search for 
X(3872) in the initial state radiation (ISR) process e + e~ — ► X(3872) r yisR 
via the decay to J/ifnr + 7r~ has yielded a null result [ID]. This strongly dis- 
favors a J PC = 1 assignment to X. On the other hand, strong evidence 
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Figure 2: The J/ip7T + n invariant mass distribution in the signal region for 
(a) B~ -> X(3872)K~ and (b) B° -> X(3872)Kg decays from BABAR. 



of the radiative decay X(3872) -> J/V7 by both BABARpl] and BellepS] 
experiments implies C = +1. This when added together with angular anal- 
ysis results from Belle [13] and CDF[H], a J PC value of 1++ for the X(3872) 
seems most likely. The only known charmonium candidate that can come 
close to this pattern of measurements is Xci> but it has an expected mass of 
3950 MeV. Clearly an alternative approach beyond the traditional charmo- 
nium model is needed to explain this state. The small width of the X(3872) 
and its proximity to the D°D*° threshold, 3871.8 ± 0.5 MeV, have led to 
the current two most-popular proposals; i.e. it could be a weakly bound 
DD* molecule-like state[15] or a diquark-antidiquark state[16]. To investi- 
gate these models, both BABAR and Belle have updated their earlier studies 
with higher statistics, and by looking into other decay modes. 

Using a dataset comprising about 232 million BB pairs, BABAR has stud- 
ied the decays B~ -> J/ij; 7r+7r'K~ and B° -> J/ip7r + 7T~K^. Figured 
shows the <//■?/> 7r + 7r~ invariant mass distributions for these two B decay 
modes [17] which are fitted to a nonrelativistic Breit-Wigner shape for the 
signal, and to a linear function for the nonresonant and combinatorial back- 
grounds. For the charged B mode, we obtain 61.2 ± 15.3 while for the B° 
mode only 8.3 ± 4.5 signal events. These yields are then translated to the 
respective branching fractions by taking efficiency corrections into account: 
Br = B(B~ -»■ X(3872)X",X -> J/^n+n-) = (10.1 ± 2.5 ± 1.0) x 10~ 6 
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Figure 3: The projections of (a) Q-value and (b) AE distribution in the 
signal region from Belle. Dots are the data points, the hatched histogram 
corresponds to the combinatorial background; the dashed histogram denotes 
the total background while the solid line is the combined fit result. 

and B° = B(B° -> X(3872)K° S , X -> J/ij, tt+tt-) = (5.1 ± 2.8 ± 0.7) x 10~ 6 , 
where uncertainties are statistical and systematic, respectively. From these 
we derive a ratio of the branching fractions, 1Z = B°/B~ = 0.50 ±0.30 ±0.05. 
We also measure the mass difference of the X(3872) state from the charged 
and neutral B decay modes, Am, to be 2.7 ± 1.3 ± 0.2 MeV. The diquark- 
antidiquark model predicts 71 = 1 and Am to be 7 ± 2 MeV. The expected 
ratio of branching fractions is consistent with our measurement, 0.13 < TZ < 
1.10 at 90 % CL, and the observed Am is consistent with zero, and with the 
model prediction within two standard deviations (er). This result seems to 
slightly disfavor the molecule model that predicts 1Z to be at most 10 %. 

Belle has reported the observation of a near-threshold enhancement in the 
D°D°7T system from both neutral and charged decays of B — > D°D°tc K [T8j. 
To determine the exact peak position as well as the branching fraction, the 
two-dimensional distribution of AE i.e. the difference of energy of the B 
candidate and the beam energy vs. Q-value (= m fl0 co T o — 2m D o — m^o) 
is fitted. Projections onto Q-value for |Ai?| < 25 MeV and AE for 6 MeV 
< Q- value < 14 MeV are shown in Figure [3] along with the total fit result. 
The observed D°D°ii mass is found out to be 3875.2 ± 0.7±?;e ± 0.8 MeV, 
where the first error is statistical, the asymmetric one is due to uncertainty 
on the 7T° energy scale and the last one is due to the uncertainty associated 
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with the D° mass. Assuming the threshold peak is due to the production 
and subsequent decay of X(3872) to D°D 7r° final state, this measurement 
is 2 a higher than the world-average value of X(3872) mass reported ear- 
lier. Corresponding branching fractions for the neutral and charged modes 
are B(B° -> X(3872)Xg,X ->_D°_D 7r°) = (1.66 ± 0.70lgp) x 10" 4 and 
B{B~ -> X{3872)K~,X -»■ D°D°7r°) = (1.02 ± 0.3lj£i) x 10~ 4 . This can 
be translated to TZ (ratio of branching fractions) to be little over one which is 
in reasonable agreement with the diquark-antidiquark model. However, we 
need more data to discriminate convincingly between these two models. 

3.2 X(3940) in double charmonium 

At the B factories exclusive B decays are not the only source of charmonium 
states. BABAR has also studied double charmonium production in the process 
e + e~ — > 7* — > J /if) cc using 124fb _1 of data[19]. Only cc states with even C- 
parity are expected in this reaction, although if there is a contribution from 
e + e~ — > 7*7* — > J/if)cc, odd C-parity states could also be produced. In this 
study, we first reconstruct a J /if) candidate via the e + e~ or decay mode 

and then calculate the mass of the system recoiling against it. Three even C- 
parity charmonium states, rj c , Xco, and rj c (2S) are observed, while there is no 
evidence for any odd C-parity state such as the J /if). The distribution is fitted 
to obtain the yield for each state, from which the production cross section is 
calculated. Due to the requirement of at least five charged tracks in the event 
for background suppression purposes, we report the product of the production 
cross section and the branching fraction to states with more than two tracks. 
The results are 17.6 ± 2.8j£f fb, 10.3 ± 2.5j£f fb and 16.4 ± 3.7l§;gfb for 
Vc, Xco and T] C (2S), respectively. The Belle measurement j2(]] in these three 
charmonium states agrees with BABAR and both are an order of magnitude 
higher than those predicted by nonrelativistic QCD[21J. However, recent 
works incorporating charm quark dynamics [22] seem to narrow down the 
discrepancy. 

An interesting addendum to these measurements is the observation of 
a new resonance by Belle [2D] at a Breit- Wigner mass of 3943 ± 6 ± 6 Me V 
and width 15.4 ± 10.1 MeV as shown in Figure HI The new state X(3940) is 
found to decay to DD*, but not to DD or ujJ/ip. The dominance of the DD* 
decay mode hints at the possibility that the A (3940) is the 2 3 Pi(cc) x'd state. 
It is natural to try the 2P(cc) since 2 3 Pj states are predicted to lie in the 
mass window of [3920, 3980] MeV with expected widths 20 - 130 MeV. The 
problem with this interpretation is, however, there is no clear-cut evidence for 
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Figure 4: The X(3940) signal in the J /if) recoil mass spectrum of the double 
charmonium production process from Belle. Black dots are data, the dashed 
curve is the background contribution while total fit result is denoted by the 
solid line. 

the 2 3 Pi(cc) Xd in the plot. Thus one is bound to suspect that the x'd should 
not be a strong signal in this reaction. This has led to alternative speculations 
that X(3940) could be the radially-excited j)". Unfortunately this also has 
its own problem as the expected rj" mass is 4040 — 4060 MeV, approximately 
100 MeV higher than the measured value. An angular analysis of the decay 
products will cast some light on the quantum numbers of X(3940) and hence 
help on unambiguously identifying the state. 

3.3 The y(3940) - is it a twin? 

Belle claims the discovery of a second resonance at 3940 MeV [23], this time 
as a threshold enhancement in the u)J/if> subsystem of the process B — ► 
KuJ/if), uj — > 7r + 7T~7r°. The reported mass and width are 3943 ± ll(stat) ± 
13(syst) MeV and 87 ± 22(stat) ± 26(syst) MeV, respectively. The new state 
is referred as V(3940) and has not been seen in the decay mode Y — > DD 
or DD* . The mass and width being similar to those of V(3940) suggests Y 
could be a radially excited P-wave charmonium. However, the uJ/if> decay 
mode is peculiar in a sense that Belle measures B(B — > KY (3940)) 13 (Y — > 
coJ/iP) = (7.1 ± 1.3 ± 3.1) x 10" 5 while we expect B[B -> K X ' cJ ) < B{B -> 
K Xc j){= (4 ± 1) x 10~ 4 }. So if we have to identify V(3940) with X ' c j, it 
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Figure 5: The J/ifm + 7t invariant mass spectrum in the range 3.8— 5.0 GeV 
and (inset) over a wider range that includes the ip{2S) state from BABAR. 



would imply B{Y — > uJ/ip) > 12%, which is quite unusual for a canonical 
cc state above the open charm threshold. 

Thus the V(3940) is something of an enigma, driving the claim of the 
Belle collaboration that it is a charmonium hybrid. The large width coupled 
with the unusual decay mode are consistent with the claim that a hybrid 
should have strongly suppressed DD, DD* or D*D* decay modes. However, 
a mass less than 4000 MeV is in conflict with the lattice QCD computations of 
low-lying hybrid spectrum. Thus more data are required before any concrete 
statement can be made on the exotic nature of F(3940). 



3.4 Observation of the y(4260) 

ISR events produced in the T(4S) energy region at the B factories serve to 
probe interesting physics occurring at lower center-of-mass energies. Moti- 
vated by this, BABAR has investigated the process e + e~ — ► J /ip n + it~^ IS R 
across the charmonium mass range, using a data sample of 233 fb -1 inte- 
grated luminosity [24] . These events are characterized by two pions, two 
leptons (electron or muon) making a J/tp candidate and a small recoil mass 
squared against the J/ip7i + ir~ system. Figure [5] shows the J/ipn + 7i~ invari- 
ant mass spectrum for the selected candidates. A clear enhancement near 
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4.2 GeV is observed in addition to the expected ip{2S) peak. No other struc- 
ture is evident in the spectrum, including the X(3872). Using a maximum 
likelihood fit in which the signal is described by a single relativistic Breit- 
Wigner lineshape, we obtain a yield of 125 ± 23 events with a statistical 
significance of 8cr (the signal is referred to as the F(4260)). Assuming the 
peak is due to a single resonance, the mass and width are determined to be 
4259 ± 8ig MeV and 88 ± 23±4 MeV, respectively. We also calculate a value 
of r(F(4260) -> e + e~) • B(Y -> J/^tt+tt-) = 5.5 ± 1.0^;?eV. Detection in 
the ISR process indicates J PC = 1 for V(4260). We have weak evidence 
for y(4260) in the exclusive B decays[TT]. In addition, it has been confirmed 
by the CLEOJ22] and Belle collaborations although at a slightly higher 
mass. 

There are number of competing theoretical interpretations floating around 
the market to explain V(4260), such as charmonium hybrid [27], hadronic 
molecule [28] etc.. However, we would need a multitude of measurements 
with higher statistics to further pin down its properties, and hence the ac- 
companying model. 

3.5 Hot of the press - the y(4320)? 

In an effort to clarify the nature of the Y (4260), BABAR has carried out a study 
of the process e + e~ — > ip{2S) 7r + 7r _ , ip(2S) — > J/ipir + ir~ using ISR events 
over a wide kinematic range from the threshold up to 8GeV(29]. The data 
represent an integrated luminosity of 272 ftT 1 recorded at y/s = 10.58 GeV, 
near the T(AS) resonance, and 26 fb -1 recorded near 10.54 GeV. First a J/tp 
candidate is reconstructed through the decay J/ip — > £ + £~ (£ = e or /i) which 
is combined with a pair of pion tracks to form ip(2S). Another pair of pions 
is then combined with the ip(2S) candidate and the resultant invariant mass 
distribution is studied. The detection of the accompanying ISR photon is 
not explicitly enforced as it is preferentially produced along the beampipe. 
Later, however, by studying the missing momentum direction we confirm the 
ISR characteristics of the event. 

Figure [6] shows the invariant mass distribution of the 2(7r + 7r~)J/ip system 
up to 5.7 GeV in the final data sample. A structure around 4.32 GeV is 
observed in the mass spectrum. To better understand this structure, we fit 
the full mass region under several hypotheses assuming the peak is due to: 
(a) V(4260), (b) '0(4415) or (c) a new resonance whose mass and width kept 
as free parameters in the fit. The probabilities of fit x 2 are determined to be 
6.5 x 10~ 3 , 1.2 x 10 -13 and 29 %, respectively. The new resonance hypothesis 
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Figure 6: The 2(ti + 7i^)J/i(j invariant mass spectrum up to 5.7 GeV for the fi- 
nal sample from BABAR. Red dots are data points while the shaded histogram 
and the curves are explained in the inset text. 

that best describes the structure has a mass of 4324 ± 24 MeV and a width of 
172 ±33 MeV, where the errors are statistical only. Given the close proximity 
to V(4260) and statistical significance, more data would clarify the nature of 
this structure. 

3.6 The Z(3930) - the last alphabet? 

This state is observed by the Belle collaboration in the two-photon mediated 
process, e + e" -> 7*7* -> DD [30] with a mass 3929 ± 5 ± 2 MeV and a 
width 29 ± 10 ± 2 MeV where the quoted uncertainties are statistical and 
systematic, respectively. The statistical significance of this claim is 5.5 o. The 
DD helicity distribution is consistent with the J = 2 expectation. Similar 
to V(3940) and V(3940) the Z seems an obvious candidate for x' C 2- The 
predicted mass and width of x' c2 are 3972 MeV and 80 MeV, respectively. 
However, if we set the x' C 2 mass to the measured value of 3929 MeV, its 
predicted width drops down to 47 MeV (reasonably close to the measurement) 
due to phase-space restriction. So at this stage we have no reason not to 
believe that the Z(3930) is the previously unseen x' c2 charmonium state. 
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4 Conclusions 



The last few years have been very exciting for haclron spectroscopy studies 
at the B factories. Both BABAR and Belle are pioneering several sensitive 
searches for new charmonium states, some of the most recent ones are sum- 
marized in this report. As it seems Xs, Ys, and Zs are seriously challenging 
our current understanding of QCD in the nonperturbative domain. Least 
of the all, X(3872) and F(4260) have made it clear that the simple quark 
model won't work for them and we need to look beyond the known horizon! 
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